One-dimensional spin-anisotropic kinetic Ising model subject to quenched disorder 
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Large-scale Monte Carlo simulations are used to explore the effect of quenched disorder on one 
dimensional, non-equilibrium kinetic Ising models with locally broken spin symmetry, at zero tem- 
perature (the symmetry is broken through spin-flip rates that differ for '+' and ' — ' spins). The 
model is found to exhibit a continuous phase transition to an absorbing state. The associated crit- 
ical behavior is studied at zero branching rate of kinks, through analysis spreading of '+' and ' — ' 
spins and, of the kink density. Impurities exert a strong effect on the critical behavior only for a 
particular choice of parameters, corresponding to the strongly spin-anisotropic kinetic Ising model 
introduced by Majumdar et al. Typically, disorder effects become evident for impurity strengths 
such that diffusion is nearly blocked. In this regime, the critical behavior is similar to that arising, 
for example, in the one-dimensional diluted contact process, with Griffiths-like behavior for the kink 
density. We find variable cluster exponents, which obey a hyperscaling relation, and are similar to 
those reported by Cafiero et al. We also show that the isotropic two-component AB —* model 
is insensitive to reaction-disorder, and that only logarithmic corrections arise, induced by strong 
disorder in the diffusion rate. 
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PACS numbers: 05.70.Ln 05.70.Fh 05.70.Jk 82.20.Wt 



I. INTRODUCTION 

The study of non-equilibrium model systems has at- 
tracted great attention in recent years. A variety of phase 
transitions have been found characterized by critical ex- 
ponents, both static and dynamic. Of special interest 
are transitions from a fluctuating active state into an ab- 
sorbing one. A wide range of models with transitions into 
absorbing states was found to belong to the directed per- 
colation (DP) universality class [HHf. Another universal- 
ity class of interest is the so-called parity conserving (PC) 
class [!, 0, IE HI • The most extensively studied model in 
this class is branching annihilating random walks with 
an even number of offsprings (BARWe) in one dimension 
(Id). The first examples of models exhibiting a PC-type 
transition were two one-dimensional cellular automata 
studied by Grassberger Q. The prototype spm-model 
for PC-type phase transitions, the non-equilibrium ki- 
netic Ising model (NEKIM) was proposed by one of the 
authors [8]. In this model domain walls (kinks) between 
unlike spins follow BARWe dynamics. 

One dimensional models showing DP transition, such 
as the contact process, exhibit pronounced changes in be- 
havior under the effect of quenched spatial disorder, even 
for small impurity concentrations (see [9( and references 
there). On the other hand, as a recent computer simula- 
tion study (lfj| has shown with great accuracy, the parity- 
conserving class appears to be highly resistant to impu- 
rities. The same study reported similar, or even stricter, 
negative results for the one-dimensional Glauber-Ising 
model 

Another aspect of t he p roblem, however, was identified 
by Majumdar et al. [12j . who introduced a specific, in- 
herent spin anisotropy (kinetic disorder) in the Glauber- 
Ising model (Majumdar Dean Grassberger (MDG) model 
in the following), These authors found, both analytically 



and numerically, a slow logarithmic factor in the decay 
of the density of kinks p(t) for t — > oo. 

The aim of the present study is to show that in certain 
Ising-like systems (possessing two absorbing states) , local 
kinetic disorder may effectively remove one of the absorb- 
ing states. The resulting, single-absorbing state system is 
sensitive to spatial disorder, similar to models in the di- 
rected percolation class. To this end we start from a gen- 
eralization of NEKIM, the nonequilibrium kinetic Ising 
model with anisotropy (NEKIMA) [HOSl, with variarj le 
kinetic disorder, and add uncorrelated spatial impurities 
to the system. We recall that without spatial impurities 
and in the zero-branching limit considered throughout 
this paper, the exponent a, which governs the time de- 
pendence of the kink density, p(t) ~ t~ a , is invariably 
a = 1/2 (apart from possible logarithmic corrections). 
The characteristic cluster-exponents, which are more sus- 
ceptible to kinetic disorder than those governing the kink 
density, differ however from their Glauber-Ising values. 
This phase will be called spin-anisotropic Glauber-Ising 
(SAGI). For definitions and details of the models referred 
to here, and the acronyms used, see Table |TJ 

Our numerical studies show that the scaling behavior 
of the SAGI phase is resistant to quenched spatial impu- 
rities, of strength po, in the small-impurity regime. Close 
to the incipient freezing of a diffusion channel, which oc- 
curs at about pp = 0.5, and which appears to be a dirty 
critical point [9(, the exponent a jumps from 0.5 to 0.21. 
On increasing po further, a Griffiths-like [l5| phase is 
encountered, with a(po) increasing with pq. The clus- 
ter properties: the mean population size n(t) ~ t v , the 
mean-square distance of spreading of spins R 2 (t) ~ t z , 
and the survival probability, P(t) ~ t~ s are also found 
to change from their SAGI expressions (Tablc|nJ) to those 
characteristic of a Griffiths phase [l5| • We find close simi- 
larities to the results on scaling in the impure contact pro- 
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TABLE I: Summary of models. Here w denotes the spin transition probabilities, A = T(l + 8), B = T/2(l — 8) 
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cess [T6| in (almost) the whole plane of the correspond- 
ing phase diagram, with the hyperscaling law remaining 
valid. Here this behavior is found in an annihilating ran- 
dom walk- type model (ARW) i.e., without kink produc- 
tion, while production takes place on the level of spins, 
in one of the spin-channels ('+' spins in the following), 
as explained in Sect. II. 

Moreover, in the present case the SAGI phase plays 
the role of a supercritical phase: active-like concerning 
'+'- cluster-behavior while frozen for the '— ' spin phase, 
and of annihilating-random-walk type for the kink den- 
sity p(t). The dual nature of global and cluster properties 
is not contradictory and can be traced back to the fact 
that the model is highly spin-anisotropic. 

The NEKIMA model coincides, at some specific pa- 
rameter values (MDG case, FigQ} with the model of Ma- 
jumdar et al. [12] mentioned above. The joint effect of 
kinetic and spatial impurities is even more pronounced at 
and in the vicinity of this point (see the phase diagram 
depicted in Fig[21 where the most spin-anisotropic point, 
the MDG-point is at the origin). In this region of the 
(p+,Po) plane the impure scaling behavior is different: 
by increasing po we have found a change in the behavior 
of the global quantity p(t): a switches from a value of 1/2 
(apart from the very late time logarithmic correction) to 
a = 1. The cluster behavior is different for the '+' and 
' — ' spins. For the '+' it appears to follow the above men- 
tioned scheme (l6j of the impure contact process: upon 
increasing po a strong tendency towards the absorbing 
'+' phase can be seen with the hyperscaling law [TtJ sat- 
isfied. For the '— ' cluster we find ?y = z/2 = 1, and S = 0, 
with logarithmic corrections in time. 

In summary, we find that by increasing spin anisotropy 
the sensitivity of critical behavior to quenched disor- 
der is enhanced. (NEKIM -» SAGI -> MDG). It is 
worth noting that in order to maintain the symmetry 
of the NEKIMA model, the disorder affecting diffusion is 
anisotropic too, therefore it is different from the problem 
considered in [1 01 ] . where for very strong disorder diffu- 
sion can be blocked completely. 

In the Appendix a further, isotropic ARW model, 
AB —¥ 0, |l[ is investigated via numerical simulations, 



under quenched disorder. In this model, which can be 
mapped onto the pure MDG model Q the effect of im- 
purities on the annihilation and diffusion rates are in- 
vestigated separately: the density decay is unaffected in 
case of the former, while logarithmic corrections arise for 
impure diffusion rates for very long times. 

The paper is organized as follows. In Sect. II. we 
define the models and the quantities to be investigated. 
Some previous results for the exponents of spreading from 
a localized source are recalled and commented on. In a 
subsection - to enlighten the different roles played by the 
'+' and '— ' spins - our model is mapped onto a reaction- 
diffusion system. In Sect. III. the effect of spatial impuri- 
ties is investigated for two different choices of the value 
of the anisotropy-parameter: the SAGI case (A) and the 
more spin-anisotropic MDG case (B). The results are dis- 
cussed in Sect. IV., while in the Appendix a further ARW 
model, the AB — s- is investigated under the influence 
of quenched disorder. 



II. SPIN-ANISOTROPIC GLAUBER-ISING 
MODEL. 

In the one-dimensional Glauber model at temperature 
zero, the most general form for the flipping rate of spin 
is ( 8i = ±1): 

r l 

w(si,Si-i,s i+ i) = ~(l+5si-iSi + i)[l- -Si(si-i + s i+ i)\ 

(1) 

Usually the Glauber model is understood as the special 
case S = 0, T — 1 and we will use these parameter values 
in the following. Processes involving the reaction 1 kink 
— > 3 kinks are introduced via the exchange rate 

W ex (s l ,S l+1 ) = ^p(l - SiS i+1 ) (2) 

(a kink is a ' — h\ or 'H — ' configuration: domain bound- 
ary between two oppositely magnetized regions). For 
negative values of S this model (called NEKIM), shows 
a line of PC-transitions in the (p ex , <5)-plane [8]. The 
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NEKIMA model is an extension with local symmetry 
breaking in the flipping rates of the '+' and '— ' spins, 
as follows. Concerning the annihilation rates (of a spin 
in the neighborhood of oppositely oriented spins) the pre- 
scription in [12[ is followed: 



'(+;—) = l, w(-,++) = o. 



(3) 



Further spin symmetry breaking is introduced in the dif- 
fusion part of the Glauber transition rate as follows. The 
transition rates 

p = w(-,+-)=w(-,-+) = T/2(l-S) (4) 

are taken as in Eq.([T]), while u>(+;H — ) and w(+; — h) 
may take smaller values: 



P+ = w(+; H — ) = w(+; — h) < p. 



(5) 



In this way, by locally favoring + spins, the effect of 
the other dynamically-induced field arising from the rates 
(EqJ3J) is counterbalanced. The spin-exchange part of the 
NEKIM model remains unchanged, Eq.([2]). 

Spreading from a localized source at criticality is usu- 
ally described by the following three quantities 



n(t) 



R 2 (t) ~ t z 



(6) 



where n(t) denotes the mean population size, R 2 (t) is the 
mean square spreading of particles (here spins) about the 
origin and P(t) is the survival probability. In most cases 
these quantities are defined for particles, in the present 
case, however, they will be used for spins. 
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FIG. 1: Phase diagram of the impure NEKIMA model (p ex = 
0, ARW case) Different regions of behavior are marked (for 
details see text) The vicinity of p+ = is not shown as not 
discussed here. 



In NEKIMA '+' and '— ' spins are not symmetric, 
therefore we have investigated two kinds of clusters. 
Namely, the development of the '— ' cluster-seed was 
started from a wholly '+' environment while that of the 
'+' cluster from a sea of '— ' spins. We will call them '— ' 
cluster and '+' cluster, respectively. The cluster critical 
exponents for the NEKIMA were reported in [3| • Their 
values at a typical SAGI point (p+ = 0.3, p = 0.5) and 
at the MDG point (p+ = p = 0.5), for the ARW-case 
discussed here, are summarized on Table HU It should be 
noted here, however, that the precision of the simulations 
in the above cited paper does not exclude the possibility 
of the presence of some log(i) corrections. 

On the level of '+' - spins the SAGI phase (1/2 = 
p > p + ) is an active one as shown by the '-(-'-cluster 
exponents on Tablellll For producing '+' spins the rate in 
Eq.((4|) with p = 1/2 is responsible, while p + annihilates 
them. Crudely speaking the probability of the '+' spin- 
generating transition process, p — p+, can be thought of 
as the quantity corresponding to the particle production 
rate in the contact process. 

As one can see due to the anisotropics the ' — ' and '+' 
spins follow different dynamics. To understand the pe- 
culiar behavior of the species we transform their motion 
into the language of reaction-diffusion models,making a 
connection with the basic models and classes [2]. 



A. Particle picture of the spin-processes 

Let X be a particle and its absence. For the ' — ' spins 
as X-es and + spins as 0-s Eq.Q gives 

X%X -> XXX (7) 

creation, while Eq.([5]) leads to the creation 

XU^XX$, 00A^0AA (8) 

From the above three, the net process remains as 

%X -> XX (9) 

Similarly Eq.Q means the coagulation processes 

0AA -> 00A, XX% -> A00 . (10) 

From eqs.© and (fTOf the generic reactions 

X <-> 2X (11) 

follow. This corresponds to the recently studied re- 
versible model, [H, EH, H3, HH, where transition to an 
absorbing phase can occur only for zero branching rate, 
which is excluded in our case, hence we don't expect an 
active-absorbing phase transition of the '— ' spins. 

On the other hand for the '+' spins as X's we get from 
Eq.© the 
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TABLE II: Cluster critical exponents at the SAGI transition point and at the MDG point. 
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spontaneous decay and from Eq.([5]) the coagulation 

0A"A" -> 00AT, XX® -> X%% , (13) 

while Eq.([4]) gives the creation by a neighbor 

XM^XXQ), mX^QlXX (14) 

These give again XX <-> X but according to Eq.(|12p here 
the spontaneous annihilation of X, X — > is also possible 
due to Eq.JSJ), which is a necessary condition for a DP 
class transition. Hence for '+' spins a phase transition 
(at finite reaction rates) is not excluded and as we will 
see later it emerges as the effect of the disorder. The 
corresponding critical behavior is similar to that of the 
disordered DP (albeit an anisotropic one). 

III. EFFECTS OF QUENCHED IMPURITIES 

At each site i the diffusion rates of kinks p = 
io(— ; H — ) = io(— ; — h) in Eq.(f5]) are modified by adding 
quenched, uncorrelated impurities with uniform proba- 
bility distribution of the form 

Pi=p + P0 *{2e(i)-1) , (15) 

where e(i) is randomly distributed in the interval < 
e(i) < 1 and < P o < 0.5. On the other hand the rates, 

which fix the (a)symmetry of the model w(+; ) = 1, 

w(— ; ++) = and p + of Eq. (Ti"5")) are kept uninfluenced. 
In this way kinetic impurities and quenched local impu- 
rities are let to act separately, though simultaneously, 
resulting in a rich phase diagram. 

Two characteristic impure behaviors are exemplified 
by the following two choices of the fixed parameters: (1) 
p + = 0.3 (to be called the SAGI-case) and (2) p + = 0.5 
(to be called the MDG-case). 

The phase diagram of the impure NEKIMA is shown 
schematically on Fig[TJ At the MDG point the '— ' phase 
gets dominating (Section III.B.). This absorbing '— ' - 
phase is separated from the SAGI phase by the MDG- 
line, a line of critical points characterized by the cluster 
exponents given in Table [TTI and where a — 1/2. This line 
is depicted here only tentatively (especially for higher po 
values), as its neighborhood is not the subject of our 
further detailed study. The diffusion-block transition 
line corresponds to the limiting value: po — 1/2, where 
the only non-fixed spin-flip probability can become zero, 
hence diffusion blocking starts. The SAGI phase in the 



vicinity of this point will be discussed in the next Section. 
It is worth emphasizing again that in the present model 
only one transition probability gets impure (one spin-flip 
rate) in contrast to a previous investigation of NEKIMCA 
with impurities 10], where all four transition rates (spin- 
flip+annihilation) have been complemented with a term 
like that in Eq. ([T5l) . We have also looked into the effect 
of impurities in one of the annihilation rates of the two 
fixing the model's anisotropy: no effect shows up in the 
vicinity of po = 0.5 in this case as can be expected. 



A. The SAGI case (p+ = 0.3) 

We will investigate the effect of impurities below the 
MDG point i.e. for p + < 0.5 at a typical parameter 
value p+ = 0.3 (see FigfT]) by simulating a system of size 
L = 3 x 10 with periodic boundary conditions up to 
t max = 10 6 Monte Carlo steps (MCS) (throughout the 
paper time is measured by MCS). For small impurity 
rates po the behavior is found unperturbed. Correspond- 
ingly, the critical exponent of the density of kinks as a 
function of time is a = 1/2 and the exponents of the '+' 
and '-' clusters are as seen on Table |TTJ By increasing 
Po, the rate of disorder, however, the spin flip rate pi in 
Eq. (|15J) can become zero in the vicinity of po = 0.5 and 
above. To interpret the results in this region it is impor- 
tant to recall [lj| that the critical behavior of the kink 
density and that of the '+' and '— ' clusters are not con- 
nected in the present model (thus no scaling law connects 
the critical dynamical exponent Z and the cluster size ex- 
ponent z). While there is no particle production on the 
level of kinks ( ARW-process) , on the level of spins there 
is production as mentioned in the previous section. The 

spin-transition rate p, for the process ( h) — > ( — h +) 

leads to the increase of + spins. This +- spin produc- 
tion is maximal at p$ = 0. The SAGI phase is an active 
one, more precisely a '+'-active one (the ' — '-spin phase is 
frozen) . At the same time in the usual sense of Glauber- 
Ising-type models, on the level of kinks, there is no active 
phase as the kink production p ex in Eq. ^ is zero. 

By increasing p$ the system stays in this '-{-'-active 
phase, neither the kink-decay exponent a = 1/2 nor 
the cluster exponents of Table [Til change. Approaching 
Po = 1/2, however, the impure spin-flip termpj, Eq. tfTS"]) . 
can become zero at certain sites and the process diagram 
becomes as shown on Fig. [5] The presence of active zones 
between otherwise frozen-in strips is clearly seen there. 
The similarity to the time-evolution pictures (Fig. 3 in 
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FIG. 2: Space-time picture at p+ = 0.3, po = 0.5 (disordered 
SAGI) starting from a random initial distribution of '+' and 
' — ' spins. The '+'-es are black 



[221] ) showing growth of DP clusters in the impure active 
and glassy phases is evident. These can be thought of as 
the much-cited active zones of the impure contact process 
0; EH HH HII inside of the otherwise inactive substrate. 
Thus we can expect similar impure critical behavior to 
take place. This means the generic presence of scale in- 
variance, and the existence of a sub- linear growth regime. 
This regime is entered now from the active one by increas- 
ing p , through the 'dirty' critical point at po = 0.5 (no 
clean critical point exists in our case, of course). In the 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 




0.0001 



1/t 



0.0002 



FIG. 3: Local exponent values of the kink density in the 
disordered SAGI phase for different values of the impurity 
strength p = 0.56,0.501,0.503,0.499,0.47,0.4 and 0.1 (top 
to bottom). 



FIG. 4: The "dirty critical point" of the intermediate SAGI 
phase entered by increasing po at po = 0.5 is demonstrated 
by this fit. The decreasing curve corresponds to the kink den- 
sity; the higher increasing curve shows the spreading radius 
R 2 , while the lower increasing one the spin number N of '+' 
spin domains. The dashed lines on top the solid ones show 
activated scaling fit of the form (JTTJl • 



vicinity of po = 0.5 the local exponent values defined as 

\n[p(t)/p(t/m)} 



*eff 



(*) = 



ln(m) 



(16) 



(where we used m = 4) behave as on Fig[3] The fixed 
point value of the density of kinks at po = 0.5 can be 
fitted with a power of ln(t) 



p{t) ~ (lnt)~ 



(17) 



with a = 1.45943 as shown on Fig. [4] This exponent 
differs from that of the known value of the infinite ran- 
domness fixed point of the 1+ld DP (a DP = 0.38197) 
[24| , not surprisingly kinks follow an anisotropic dynam- 
ics here. The excellent fit with Eq. (fl~7|) suggests that 
this point is a dirty critical point, which is reached from 
the 'active' phase (here from the '+ '-active phase). The 
value of a is a new critical exponent. For higher values 
of impurity strength a Griffiths-like phase with 



P (t) ~ r 



*(j>o) 



(18) 



is entered. Though in this region the present formulation 
of the impure problem already goes beyond its range of 
validity, however the tendency of the impure behavior is 
well exhibited in our simulations (e.g. a(pa) increases 
steadily as expected in a Griffiths phase, etc.). Never- 
theless, we refrain from giving more details. 

The time dependences of the '+' and '-' cluster quan- 
tities (Eqs. also show a behavior characteristic of an 
impure phase transition. The behavior of the '-' clus- 
ter shows an effect of 'melting': r\ and z/2 move from 
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to some finite values, upon increasing po, while S = re- 
mains the same. The '+' cluster behavior however, shows 
a scenario of a transition to an absorbing state, similarly 
to the contact process investigated by Cafiero et al. [16| • 
However the direction is opposite: in for small values 
of disorder (called p in their paper) an absorbing phase 
is present, while the active is reached at p = 1. 

We plot on Fig. [5] the '+' cluster behavior, namely the 
local exponent values of the survival probability defined 



5eff(t) 



-\n[P(t)/P(t/m)} 
ln(m) 



(19) 



for different values of disorder. One can clearly see a 
turn from the SAGI-active behavior (S = 0) towards a 
Griffiths-phase-like non-universal cluster scaling by in- 
creasing S in the vicinity of po = 0.5. 

At po = 0.5 (dirty critical point of the '+' spins) an 
ultra-slow, activated scaling of the form Eq. (fPTj) can be 
observed in the cluster data (see Fig. 2]). Fitting with 
the activated scaling forms 



d(t) ~ (lnty b ,R 2 (t) ~ (lnt) z ,N(t) ~ (Int) 7 



(20) 



for the '+' data resulted in: 5 = 0.060(6), z = 3.88(1) 
and f\ — 3.37(2). These values again differ from those of 
the strong disorder fixed point of DP, probably due to 
the anisotropies of the processes. 

Our simulation data for the mean size of the '+'-cluster 
R(t) depicted on a double-logarithmic scale is shown on 
Figj6] The linear SAGI behavior crosses over to decreas- 
ing values of z/2 by reaching the po — 0.5 transition 
point (sub-linear behavior, for similar behavior see Ta- 
bles I. and II. in 0) 



B. The MDG-case (p+ = 0.5) 

The MDG point can be found at the origin of the 
phase diagram (FigfT]). At these parameter values the 
impure system decays into a fully compact state, where 
all spins are flipped to '— ' but different types of kinks 
exert marginal perturbation on each other [25| . In more 
detail (l2J found analytically that the magnetization de- 
cays as 



m(t) = — 1 + const/ ln(fet) 
while the kink density behaves as 

p(t) ~t- 1/2 ln(6t) 



(21) 



(22) 



(where b is constant). For the time being - since 
much stronger (power-law-type) effects enter due to the 
quenched spatial impurities - we disregard the ln(t)-type 
of correction and show only the leading singular behavior 
of the cluster properties (see Table [TTJ) . 

By introducing impurities of strength po = 0.3, the 
space-time evolution of kinks looks as can be seen on 
FigfT] The absorbing nature of the '+' phase is apparent. 
The fully '— ' spin configuration sets in quickly (compare 
this figure with Fig. 1. of [14], where single '— ' spins 
move via random walk in the see of '+'-es.). 

Results of high-precision simulations are shown on 
FigfS] For po < 0.5 one can see a — > 1 in the long 
time limit. This behavior can be ascribed to the choice 
of the impurity, which acts against one type of diffusion 
(which favors '+' spins). As a consequence the '— '-es 
can spread better than in the impure case. The effect of 
impurities on the '— ' clusters is illustrated on Fig. [9]for 
the rj and z/2 exponents. The local slope curves of r\ and 
z/2 are almost the same, they approach the asymptotic 
value rj = z/2 = 1 from below with an upward curvature 
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FIG. 7: Space-time evolution in the disordered MDG model 
(p+ — 0.5, pa = 0.3) from the random initial configuration. 
The + (black) and — (white) spin clusters are separated by 
kinks 
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FIG. 8: Local exponent values of p(t) for different values of 
po = 0.08, 0.18, 0.25, 0.45, 0.5 (from top to bottom) indicating 
a = 1-type power-law behavior. 



(<5 = remains constant). This kind of effective exponent 
behavior usually corresponds to logarithmic correction to 
scaling. In case of the general form for N(t) 




N{t) = t"In 6 (t) 
the effective exponent behaves as 
dln(N(t)) 



Veff 



= V 



dln(t) 

Therefore a logarithmic fit of the form 
a + b/ ln(i) 



ln(t) 



(23) 
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FIG. 9: Negative cluster exponents rj and z/2 in the impure 
MDG case for po — 1, 0.56, 0.35, (in pairs, from top to bot- 
tom). 5 = (P(oo) — 1) in the whole interval. The dashed 
lines are logarithmic fitting of the form 



IIII1 The clusters are compact, of course, thus the hyper- 



TABLE III: Logarithmic fitting results for the ' — ' spins clus- 
ters at the MDG point for various disorder via Eq. (|25p . 





Po = 1 


po = 0.56 


po = 0.35 


7, = z/2 

b 


1.001(1) 
-1.511 


1.004(5) 
-1.825 


0.95(6) 
-1.87 



scaling law is satisfied. 

To summarize, instead of the weak logarithmic ten- 
dency towards the dominance of the '— '-es as for the 
pure MDG model the impurity pushes the system into 
a more anisotropic decay, characterized by the kink de- 
cay exponent a = 1. The depletion of the '+ '-cluster is 
illustrated on FigsfTUl and [HI The '-(-'-cluster dies out 
apparently in a power-law manner with the exponents 
X] = — 1, 8 = 1 and z/2 — 0. These scaling behaviors 
suggest a transition from one absorbing phase (MDG) to 
an other one as the effect of quenched impurities, which 
increases anisotropy. 

It is worth noting that by switching on the spin- 
exchange term (p ex ) in the pure case a quicker 
dominance of the '— '-es takes place for the same 
parameter values of the model as here. In that case, 
no power-law regime is visible (except for some early 
time transient region and for very low values of the 
kink production p ex ) and the absorbing phase sets in 
exponentially fast typically (see [l3|]) due to the faster 
destruction of metastable '+' domains. 



has been applied for the local exponent curves and the 
values obtained both for r\ and z/2 are shown in Table 
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0.001 0.002 

1/t 

FIG. 10: Local exponent values of the '+' cluster exponent 
z/2 in the impure MDG case for p = 0, 0.15, 0.20, 0.35, 0.45 
(from top to bottom). 
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FIG. 11: Local exponent values of the '+' cluster exponent 8 
in the impure MDG case for p = 0, 0.15, 0.20, 0.35, 0.45 (from 
top to bottom). 



IV. DISCUSSION AND CONCLUSIONS 

In this paper the problem of the effect of quenched 
impurities on one-dimensional non-equilibrium Glauber- 
Ising-type models has been investigated numerically. 
Our previous study of the same problem, without spin- 
anisotropy, led to the conclusion of a minimal of effect 
[To| . Here the emphasis has been on introducing spin- 
anisotropy, and to give insight into the the joint effect 
of kinetic and quenched spatial anisotropies. We have 
restricted our studies to the zero branching rate limit 
of the model (ARW-case). Depending on the ratio of 
anisotropies in the diffusion and annihilation channels 
(rates) a variety of behaviors appear when the quenched 



randomness is turned on, with the inference that some 
amount of kinetic anisotropy can significantly enhance 
the influence of quenched spatial anisotropies. 

This behavior is different from that found for isotropic 
NEKIM [10( , where reaction and diffusion are completely 
blocked for very strong disorder. The effects of such 
blocking are investigated in the context of the AB — ► 
model; only the diffusive ones are found to be marginally 
relevant. 

Summarizing the results for exponents 77, 5, and z (Ta- 
ble HU, for '+' spins in the case p = p + = 1/2 (MDG+) 
we have the (0, 1/2, 1) spreading fixed point exponents of 
the XX — > X diffusive coagulation process, because due 
to Eq.([3]) the XX <-*• X is unbalanced, favoring XX — » X 
annihilation of '+'-es and the proliferation of '— '-s. In 
the p + < p = 1/2 kinetic disorder case (SAGI+) we 
found [lj] a novel kind of anisotropic fixed point (1, 0, 2), 
which is the same what one obtains if quenched disor- 
der is added to (MDG-). This can be understood by 
looking closer on the effect of these disorders on the spin 
species. Both kind of disorders enhance the production 
of one species at the expense of the other. While in 
case of (SAGI+) the '— '-es are suppressed, for the (diso 
MDG-) the same thing happens with the '+'-es (due to 
Eci. (fl~5|) and the other spin domains grow linearly with 
77 = z/2 = 1. 

However due to Eq.® the handicapped species end 
up in a different way. In case of SAGI the '— '-s can sur- 
vive resulting in (0,0,0) as t — > 00, slowing down the 
domain wall decay (resulting in a — 0.5). In case of the 
diso-MDG the lonely '+'-s can die out completely, char- 
acterized by 77 = — 1 mass and 6=1 survival exponents, 
enabling the linear growth of kinks. Therefore a = 1 can 
be observed here. 
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VI. APPENDIX 

In this appendix the effect of quenched disorder in a 
further ARW type model, the AB — * one, will be inves- 
tigated. Unlike SAGI, this model is free of anisotropies. 
In our previous paper [loj we investigated the inactive 
phase of the cellular automaton version of the nonequi- 
librium kinetic Ising model (NEKIMCA) model with 
quenched disorder. In that model annihilating random 
walk of domain walls (kinks) AA — > is the dominant 
process. We have not found any effect of the weak dis- 
order on the density decay of kinks. If the disorder was 
so strong that diffusion could be blocked at certain sites 
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(the system is split into independent blocks of finite sizes) 
a new fluctuating phase with finite steady state density 
emerged as the consequence of the parity conservation 
of kinks. Along and just below the blocking transition 
line the disorder caused logarithmic slow-down correc- 
tions to the p oc t~ x / 2 impure law of the ARW process 
[2(| and continuously changing exponents (a < 0.5) for 
strong disorder. Those results are i n ag reement with an- 
alytic results for diffusion disorder [27[ and annihilation 
disorder g| of the ARW. 

0.4 | , , , , , , 1 



0.35 - 



0.25 - 




°' 2 10' 10 2 10 3 10" 10 5 10 6 10 7 10 s 

t 

FIG. 12: Density decay of the annihilation disordered AB — > 
model in case of RI conditions for ea = 1, 0.98, 0.95 (top to 
bottom) 

The anisotropic reactions of the NEKIMA model [l3| 
result in a two component reaction dynamics of kinks 
destroying the PC class critical behavior. Without kink 
production (p ex = 0) one can map onto the AB — > 0, 
however, since A and B species are domain walls of the 
spins they are arranged alternately along the Id system. 
Consequently they can always meet their nearest neigh- 
bors (no hard-core reactions take place as in [2!|) and 
one observes the usual p oc i -1 / 2 decay instead of the 
p oc i" 1 / 4 law of the two component model (30l |. 

In this appendix we show results for the quenched dis- 
ordered AB — > model to complement results for the 
quenched disorder NEKIMA model. To our knowledge 
this has not been studied before. We have investigated 
the effect of disorder in diffusion and annihilation sep- 
arately. A random sequential simulation program has 
been run in Id systems of size L — 5 x 10 5 with periodic 
boundary conditions and single particle occupation re- 
striction. We have investigated the density decay for ini- 
tial conditions with randomly (RI) and pair-wisely (PI) 
distributed ^4-s and B-s of full lattice occupancy. 



One elementary MC step is built up as follows. A 
particle and a direction is chosen randomly. If the nearest 
neighbor (nn) site in the selected direction was empty the 
particle is moved there with probability 1 — en (where 

to is the diffusion disorder strength). If the nn site was 

0.4 | , , , 1 



0.35 - 



0.25 - 




°' 2 10° 10 2 10 4 10 6 10 s 

t 

FIG. 13: Density decay of the diffusion disordered AB — > 
model in case of RI conditions . Different curves correspond 
to to = 1,0.95,0.5 (top to bottom). 



occupied with a particle of a different type we removed 
both particles with probability 1 — €a- The time (t) is 
incremented by At = 1 jn v (where n p is the total number 
of particles) and the density decay was followed up to 
t max = 10 s MCS. 

In case of try = we did not see any effect of the 
disorder in the annihilation up to < 1 both for RI 
and PI conditions (see FigfT2]for RI). The density decay 
is completely insensitive to the quenched randomness of 
the annihilation reaction rate. 

On the other hand for disorder in the diffusion (hop- 
ping) probabilities (e^ = 0, en < f) one can see log- 
arithmic corrections to the scaling of the impure case 
for en — > 1 (see Fig. [T2]for RI). Following a sharp ini- 
tial decay, when local correlations are built up the evo- 
lution tends towards the p oc <~°- 25 law, but as in case of 
the ARW model logarithmic slow-down appears for very 
late times and strong disorder (already without diffusion 
blocking). It has been tested by simulating on different 
sizes (L) that the turns in the curves are not artifacts of 
finite system sizes. Since the basic mechanism of anni- 
hilating random walk here AB — ^ is the same as for 
the AA — > model, we think that the relevancy of diffu- 
sion disorder compared to the annihilation disorder can 
be observed in the inactive phase of the NEKIMCA [l(| 
as well. 
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NEKIM=Glauber-Ising model + n.n. spin exchange (eqs.(l) and (2)) 

pex=0 pex>0 
Annihilating Random Walk Branching Annihilating Random Walk 
(ARW) *" (BAR We) 

GLAUBER-ISING CLASS PARITY CONSERVING CLASS 

( PC - transition if 5<0 ) 

+ SPIN ANISOTROPY = NEKIMA 



w(+;- -)=1, w(-;+ +)=0 (MDG rule) 
w(+;+ -)=w(+;- +)= p + 

w(-;+ -)=w(-;- +) = p=l/2 

p^Fp MDG-case 
p + < p SAGI-case 



(DP -transition) 



NEKIMCA 

cellular automaton version of NEKIM (pex=0) 

update of spins: synchronous 
leads to kink production and PC transition 



